On the Pionner Anomaly 

Kostadin Trencevski 

Faculty of Natural Sciences and Math., Ss. Cyril and Methodius Univ., P.O.Box 162, 1000 

Skopje, Macedonia 

Jq ■ kostatreOiunona . pmf . ukim . edu . mk 

O , and 

(N 

p ■ Emilija G. Celakoska 



< 



Faculty of Mechanical Engineering, Ss. Cyril and Methodius Univ., P.O.Box 464, 1000 

Skopje, Macedonia 



cemilSerebl . mf . ukim . edu . mk 



> 

00 
\D 
O ■ ABSTRACT 

o 
^ . 

O . In this paper an explanation of the Pioneer anomaly is given, using the hy- 

O . pothesis of the time dependent gravitational potential. The implications of this 

hypothesis on the planetary orbits and orbital periods are given in section 2. In 
bX). sections 3 and 4 we give a detailed explanation of the Pioneer anomaly, which 

corresponds to the observed phenomena. 

X, 

j^ . Subject headings: Pioneer anomaly, Hubble constant, gravitational potential 



1. Introduction 

In 1998 an interesting paper by Anderson et al. (1998) was published, with indication 
that there is an unexplained frequency shift in the navigation of the spacecraft Pioneer 
10 and 11. Namely, the detailed analysis of the frequencies of the round-trip signals from 
the spacecraft Pioneer 10 and Pioneer 11, yield to a slight disagreement with the expected 
frequencies. This small shift of the frequencies accumulates in time. Later a very detailed 
analysis of this phenomena was given by Anderson et al.(2002). Firstly, we give a short 
history about this effect, following the references (Anderson et al. 1998, 2002; Turyshev et 
al. 2004, 2005, 2006). 
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The spacecraft Pioneer 10 was launched on 1972 March 2, while Pioneer 11 was launched 
on 1973 April 5, and their aim was to explore the outer Solar System (Turyshev et al. 2006). 
They are the first spacecrafts that left the Solar System. The mentioned phenomena was 
detected in 1980 for the spacecraft Pioneer 10 for the first time, when it was about 20 AU far 
from the Sun, and the influence of the neighboring planets was negligible on the detected shift 
frequency. This phenomena was well examined for Pioneer 10 on the heliocentric distances 
20-70 AU. Similar results were also obtained for the spacecraft Pioneer 11, but, unfortunately, 
after 1995 November the contact with Pioneer 11 was lost, and the research with Pioneer 
11 was stopped. The last signal from Pioneer 10 was obtained in April 2002, when it was 
about 80 AU from the Sun. 

The unexplained blue shift of the frequency is denoted by ap and it is introduced by 
the formula (15) in (Anderson et al. 2002), i.e. 



r / \ / M 2cipt 
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where z/q is the reference frequency, the factor 2 is because of the two way data, Vmodei is 
the modeled velocity of the spacecraft due to the gravitational and other forces and Uots 
is frequency of the re-transmitted signal observed by DSN antennae. Since ap is in the 
same units as the acceleration, it is also called anomaly acceleration and it is interpreted 
as if there is an acceleration toward the Sun of magnitude ap ~ (8.74 ± 1.33) x 10~^cm 
s~^. The existence of such an acceleration is almost impossible, because it would change the 
planetary orbits, while the astronomical observations do not show such change (lorio 2006). 
So it would be more convenient to say "Pioneer anomaly" than "anomaly acceleration". It 
is mentioned in (Anderson et al. 1998, 2002) that without using the acceleration ap, the 
anomalous frequency shift can be interpreted as a "clock acceleration" —at = —2.8 x lO^^^s 
s-^ 

There are several attempts to explore the possible sources of the Pioneer anomaly. The 
original efforts, including classical explanations (effects from sources external to the space- 
craft, study of the on-board systematics and computational systematics) of this anomaly are 
not strong enough to explain the observed Doppler shift (Turyshev et al. 2006). So, the 
question of a possible new physics arises. There are several attempts to explore the possible 
sources of the Pioneer anomaly via non-classical Newtonian mechanics and general relativ- 
ity. The ideas are quite different and here we mention some of them: Newtonian gravity as 
a long wavelength excitation of a scalar condensate inducing electroweak symmetry break- 
ing (Consoli & Siringo 1999), time dependence of the gravitational constant G (Sidharth 
2000), higher order (or Yukawa-like) corrections to the Newtonian potential (Anderson et 
al. 1998), a scalar-tensor extension to the standard gravitational model (Clachi Novati et al. 
2000), interaction of the spacecraft with a long-range scalar field coupled to gravity (Mbelek 
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& Lachieze-Rey 1999; Cadoni 2004), a length or momentum scale-dependent cosmological 
term in the gravitational action functional (Modanese 1999), a local solar system curvature 
for light-like geodesies arising from the cosmological expansion (Rosales & Sanchez-Gomez 
1999), influence of the cosmological constant at the scale of the Solar System (Nottale 2003), 
a 5-dimensional cosmological model with a variable extra dimensional scale-factor and static 
external space (Belayev 2001), a superstrong interaction of photons or massive bodies with 
the graviton background yielding a constant acceleration proportional to the Hubble con- 
stant (Ivanov 2001), using an expanded PPN formalism (Ostvang 2002), explanation using 
the Brans-Dicke theory of gravity (Capozzielo et al. 2001), explanation using braneworld 
theories (Bertolami & Paramos 2004), and so on (Foot & Volkas 2001; Guruprasad 1999, 
2000; Wood & Moreau 2001; Crawford 1999; Ingersoll et al. 1999; Milgrom 2001; Munyaneza 
& Viollier 1999; Ranada 2005). 

In this paper, we give new arguments and we upgrade the papers (Trencevski 2005a, 
2006) concerning the hypothesis the time dependent gravitational potential. In section 2 we 
give some results about the time dependent gravitational potential and in section 3 we will 
return to the Pioneer anomaly. The mathematical calculations will be shortened, while the 
physical interpretation will be emphasized. 



2. Time dependent gravitational potential 

In the recent paper (Trencevski 2005a) and also in (Trencevski 2006), an idea about 
linear change of the gravitational potential in the universe was developed. This linear (or 
almost linear) change is probably a consequence of the change of the density of the dark 
matter in the universe, but we shall consider only the consequence of this phenomena, not 
its origin. It can be interpreted in this way: the whole universe is inside of a shielded labora- 
tory which falls freely in a gravitational field. Then, for a small time interval, the change of 
the gravitational potential can be considered to be almost linear. For simpler interpretation 
of the induced effects, we shall consider two observers: observer A in the shielded laboratory, 
where the gravitational potential changes linearly, and observer B outside the shielded lab- 
oratory, far from the gravitational field, where the gravitational potential is constant. Thus 
the proper time for the observer B is uniform, i.e. not accelerated. 

The coefficient of the linear change of the gravitational potential is uniquely determined, 
so the observer A measures the Hubble redshift: The light signal, which starts from a distant 
galaxy, comes to the observer A after a long period, when the time is much faster, and so A 
observes a redshift. 
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If a light signal starts from a star with a frequency uq, then, according to the general 
relativity, after a period t = — its frequency will be 

/ RdV\ 

It is assumed here that the sign of the gravitational potential V is such that the gravitational 
potential V is larger near the massive bodies. In comparison with the Hubble law, 

where H is the Hubble constant, H ^ 70 km s^^ Mpc^^, we obtain that 

^ = -c^H ^ -2 X 10=^cm2 s-l (2.1) 

dt ^ ' 

Notice that this linear change of the gravitational potential in the universe is the same as 
the change of the Earth's gravitational potential in a lift which moves toward the high floors, 
with velocity 2 cm s^^. 

Now, accepting this linear change given by (2.1), the simplest application is the expla- 
nation of the Hubble redshift. Thus, the main reason for the Hubble redshift is the linear 
change of the gravitational potential, while the Doppler effect has a minor role. 

Further, we shall consider the influence to the planetary orbits and orbital periods 
(Trencevski 2005a, 2006). 

Let us denote by X, F, Z, T the coordinates in the shielded laboratory observed by the 
observer B, and let us denote by x, y, z, t the coordinates observed by the same observer B 
but in his close neighborhood. This can be interpreted in the following way. Assume that in 
the shielded laboratory there is motion of a particle around a body under the gravitation. 
Then X, y, Z are functions of T on the trajectory. Assume that, simultaneously, the same 
happens out of the shielded laboratory in a close neighborhood of the observer B. Then, 
for this simultaneous trajectory, a:, y, z will be functions of t. The observer i? is in a good 
position to compare these two trajectories. 

The coordinates x^y^z^t will be called "normed coordinates", because the time is uni- 
form and the dimensions of objects in a neighborhood of B are static. According to the 
general relativity (GR), up to the first postnewtonian approximation we accept that 



dx=[l + — \ dX=\l + tHjdX, (2.2) 

dy= (l + ^) ^dY = (l+ tH)dY, (2.3) 
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dz= (l + ^) ^dZ = (l + tH^dZ, (2.4) 

dt= (l- ^) ^dT = (l - tH^dT. (2.5) 

Assuming that (2.2), (2.3), (2.4), and (2.5) are satisfied, it is not necessary to speak 
about the time dependent gravitational potential. Indeed, if we accept (2.2), (2.3), (2.4), 
and (2.5) as axioms, and neglect H^, then tH can be replaced by TH, (1 + tH)~^ = 1 —tH, 
and so on. Further on, the Hubble redshift is a direct consequence from (2.5). 

From (2.2), (2.3), (2.4), and (2.5) we obtain 

fdX dY dZ\ fdx dy dz\ 

and by differentiating this equality by T we get 

^d^X d^Y d^Z\ _ 

KdT^'lm'^m) ~ 

^ \'d^''d^''dt^) ~ \'d^''dl^''dl^) ~ y'dT'dT'dr)- ^ ' 

In normed coordinates x, y, z, t there is no acceleration caused by the time dependent gravi- 
tational potential, i.e. we can put there H = 0. Thus, according to the coordinates X, Y, Z, T 
there is an additional acceleration 

^.-^(d?x d^y dh\ ^frjdX dY dZ\ 

It is very important to notice that the right sides of (2.2), (2.3) and (2.4) are not total 
differentials, so there does not exist a functional dependence between the "coordinates" 
X, y, z, t and X, Y, Z, T. This shows that we deal with unholonomic coordinates. The required 
functional dependence between x, y, z, t and X, F, Z, T can exist only along a chosen curve, 
and it will depend on the choice of the curve. This yields to a violation of the strong 
equivalence principle, i.e. to one part of it, but on the other hand we shall see that the 
previous equations (2.2-5) yield to results compatible with the experimental measurements. 
The violation of the strong equivalence principle is explained in more details in (Trencevski 
2005b). 

In order to see the influence of the constant H to the planetary orbits, we are looking 
for a functional dependence of the form 

a; = (1 + \tH)X, y = {l + \tH)Y, z = {I + \tH)Z, dt = {1 - fj,tH)dT, (2.8) 
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A = const, and /i = const., that will lead to the same equality (2.7). As a consequence from 
(2.8) it is 

= (1 - (A + 2,)TH){^, -^.-^) - (2A + M)i/(-. -. -). (2^9) 

Comparing (2.7) and (2.9) we obtain A = | and /^ = |- Now we give the main consequences. 
In order to distinguish the influence of the time dependent gravitational potential, we assume 
Keplerian orbits for H = 0. 

1. From (2.8) it follows that the quotient ©2 : Oi of two consecutive orbital periods 
observed from B is equal to 1 + jiQH = 1 + |0-f^- This shows that each next orbit has a 
prolonged period for a factor 1 + ^QH . But according to the same observer B, the time 6 
for the observer A in the shielded laboratory is prolonged for a factor 1 + QH according to 
(2.5). Thus, according to the observer A, each next orbit is prolonged for the coefficient 

This is a consequence of the unholonomic coordinates. Formula (2.10) can be verified for 
the orbital periods of double stars (Trencevski 2006). From (2.10) we obtain 

A = ^nif, (2.11) 

where P^ is the orbital period of the pulsar in a binary system. Formula (2.11) is tested 
for the binary pulsars B1885+09 (Kaspi et al. 1994) and B1534+12 (Stairs et al. 1998, 
2002). which have very stable timings, and the results are satisfactory. Indeed, formula 
(2.11) together with the infiuence of decay of the orbital period caused by the gravitational 
radiation and a non-gravitational infiuence of kinematic nature in the galaxy, yield together 
to the measured value of Pb (Trencevski 2005a, 2006). Notice that if we neglect the infiuence 
from the equation (2.11), then the change of the orbital period caused by the gravitational 
radiation would not fit with the experiments. 

2. According to the observer B, the distance of any planet to the Sun does not change 
according to the Keplerian law, but at each moment the Keplerian distance to the Sun should 
be multiplied by the factor 1 — \HT = 1 — \HT . On the other side each distance in the 
shielded laboratory is observed by B to be shortened for the coefficient 1 — HT. Thus, at each 
moment the Keplerian distance to the Sun is observed by A as multiplied by the coefficient 
\-HT = 1 + 3^^- Here, the distances are assumed to be measured via laser signals, as for 
the LLR. 
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3. The planetary orbits are not axially symmetric and according to the observer B 
(Trencevski 2006) the angle from the perihelion to the aphelion is tt— ^^^^^— ^= tt— ''^ }^^~^ , 
while the angle from the aphelion to the perihelion is tt + ^^Qvi~g =tx + gevi-e ^j^ere 
is the orbital period and e is the eccentricity of the orbit. This follows from (2.8). Indeed, 
this happens according to the observer B who observes that the distances are shorter for 
coefficient 1 — \HT. But according to the observer A, who observes that the distances are 
increasing for the coefficient 1 + \HT, the angle from the perihelion to the aphelion is equal 
to 7r+ ^^'^Z.^'^ 1 while the angle from the aphelion to the perihelion is equal to tt— ^'^®^J~^ . 

4. Using (2.8), or the previous discussion, the time dependent gravitational potential 
has no influence on the perihelion precession. 

5. In (Trencevski 2006) are also considered the increasing of the orbital period of the 
Moon, the distance to the Moon and the change of the average Earth's angular velocity. 
These quantities depend on tidal dissipation and also on the time dependent gravitational 
potential. Including the influence of the time dependent gravitational potential, the discrep- 
ancies among the previous three changes become much smaller (Trencevski 2006). 



3. The Pioneer anomaly 

Now we shall see how the time dependent gravitational potential can explain the Pioneer 
anomaly. Firstly we raise some questions of such kind, that any viable theory that gives an 
explanation of the Pioneer anomaly should give answers. 

1. Why the considered anomaly is detected only for hyperbolic trajectories? 

2. Why the anomaly becomes observable only for heliocentric distances larger than 10 
AU? 

3. Why ap on small distances (Pioneer 11) is much smaller than its value on distances 
larger than 20 AU, where it is almost a constant? 

4. Why the induced "anomaly acceleration" toward the Sun is found only for spacecrafts, 
but not for the trajectories of planets? 

We shall find the frequency v = Vobs of the re-transmitted signal received by the DSN 
antennae. It will be calculated according to observers A and B. Let us denote by R the 
distance from the spacecraft to the Sun. 

The required frequency caused by the motion of the spacecraft, i.e. the Doppler effect. 



observed by the observer B is given by 



V 



t'o 



1 



2dT 

c 



(3.1: 



where z/q is the initial frequency. This is just the Doppler redshift for the outside motion 
of the spacecraft. Here dR is the increment observed from i?, and according to (2.2), (2.3), 
and (2.4), 

dR={\- tH)dr, (3.2) 

and dT = {l+tH)dt according to (2.5). According to the observer A (concerning the Doppler 
effect), dt should be replaced by dt{l + \tH) because of the following reason. The frequency 
Vobs in (1-1) is modeled as a function of time and there is a dilation of the orbital time for 
coefficient 1 + ^tH observed from A, which also has influence on ap. This is equivalent to 
the statement that the time dT in (3.1) is the "trajectory time", i.e. the time needed for 
the spacecraft to move from one point to a close neighboring point on the trajectory, which 
is given by 

dT = {l + -Ht\dt. (3.3) 

Now, substituting the values of dR and dT from (3.2) and (3.3) into (3.1), according to the 
observer A (3.1) takes the form 



V 



1^0 



dr 

1-2^(1 



3 



(3.4) 



Moreover, according to the observer A, in (3.4) will appear, also, an additional expression 
"^'^^ which is analogous to the redshift observed from the distant galaxies. Indeed, the time 
which passes the light signal in two directions is just — . Hence 



t'o 



f? — 

c c 



Im 

3 



(3.5) 



Note that in this expression H^ is neglected. Also, the relativistic corrections can be 
neglected here, because the unmodeled shift frequency is free of the relativistic corrections. 

Notice that the Pioneer anomaly is modeled without implementing the Hubble constant 
H and also the measurements are performed according to the observer A, for whose point 
of view the time dependent gravitational potential is applied. Instead of this formula, Uobs is 
modeled by putting H = 0, and according to (1.1), the unmodeled value ap appears. There, 
for the round-trip time period t = — of the signal, we obtain 



R 



V 



^0 



dr . 
dt 



l + 2ap--2^ 



(3.6) 
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Both formulae (3.5) and (3.6) are deduced according to the observer A. Comparing 

them, we obtain 

/7dRT \ , , 

which gives the required expression for ap. Notice that for an arbitrary initial value of the 
time T (or t), ap ^' ^cH. 

At this moment it is undetermined what is the initial value for t (or T). This problem 
will be considered in the following section. 



4. Possible values of the initial time values 

There is no global determination of the initial value of t, i.e. of T, because the physical 
laws must be independent of that choice. Still, we may consider different cases and make 
suitable choices. If the spacecraft moves on an elliptical orbit, the independence of the 
physical laws from the initial value of the time can be interpreted in the following way. 
Suppose that the orbital period of the spacecraft is 0. Then after time t + 0, from (3.5), 
we obtain 



t'o 



M .f/. 7 



l-2H--2^(l--H{t + Q) 
c c \ 3 



Now, in comparison with (3.5), we obtain that after the period when the spacecraft 

dr 
dt 



reaches the considered initial position - distance R to the Sun with the same velocity ^ at 



that position, a small blueshift will appear, such that 

14 — 

Au = —^QHuo. 

3 c 

But, a good interpretation for hyperbolic trajectories, without periodicity, is much more 
difficult to find. The interpretation in that case is given only by (3.7), but the biggest problem 
comes from the determining of Vmodei- 

Now let us calculate /?, T, and ^ as functions of one variable {R or T). It is sufficient 
to solve the problem within Newtonian mechanics. If the hyperbola of motion is given by 

2 2 

% — Ij = 1, then we choose -u as a parameter and obtain in final form 





/ 


R{u) = y 


7 / ^2 \2 , 


a ( 





(4.1) 
(4.2) 



-10- 



dR, ^ C l + #-e"' 



,u) = u , , 4.3 

di a 1/ I 7 \2 

l + t-at] +m2I 



where e is the eccentricity of the hyperbohc trajectory, C = R^^ = const ^ and K is the 



dt 

required constant. Substituting R, T and ^ from (4.1), (4.2), and (4.3) into (3.7) we obtain 
ap as function of u (and hence R) and the unknown constant K. Notice that ap does not 
depend on the constant C and thus we may put C = 1. Namely, we obtain 



ap = \-U • ^—;rz , ^9^ T~, T^ — 1 \cH+ 



l + f-ae) +u^ 



e — 



Vu'^+b'^ 



o2 -2 






3 a 



2 ^2 . ,- 



1 + ^ - ae ) +u 



e — 



cH. (4.4) 



Notice that i^ = means that measuring the time starts when the spacecraft reaches the 
perihehon position, while K > means that measuring the time starts before the spacecraft 
reaches the perihelion position. The values of K may be different for Pioneer 10 and Pioneer 
11. Moreover, K may depend also on the interval where ap is considered. Indeed, these values 
are hidden in the complicated calculations of modeling Umodei- Probably K corresponds to a 
such value, that ap seems to be almost a constant in the chosen interval. 

The values of a, b and eccentricity e for Pioneer 10 are given by (Anderson et al. 2002) 

a = 1033394633 km, b = 1463395947 km, e = 1.73359601 

and the values for Pioneer 11 are given by (Anderson et al. 2002) 

a = 1218489295 km, b = 2316289260 km, e = 2.147933251 

The values of u and R are given in AU. Notice also that if we take C = 23.2609876 
for Pionner 10 and C = 33.906372 for Pioneer 11, then the values of t are given in years. 
The values of ap of Pioneer 10 for i^ = 15 are given on figure 1, and the values of Pioneer 
11 for K = 30,45,60 are given on figure 2. Notice that for Pioneer 10 -ft' = 15 means that 
measuring the time started 235.5 days before the spacecraft was at the perihelion. Notice 
that for this value of K, ap = 1.18(1 ± 1.34%)ciJ while the spacecraft was between 20 
AU and 70 AU. Moreover, assuming that the mean value of ap is 8.74 cm s~^, then it 
corresponds to H = 76.2 km s~^ Mpc^^. This is in accordance of the observed values of ap, 
its deviations and the measured value of H, having in mind that about 10% of the observed 
value of ap may be caused by another minor reasons mentioned in section 1. For Pioneer 
11 we have a similar situation. For i^ = 45 it means that measuring the time started 485 
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days before the spacecraft Pioneer 11 was at the perihehon. For this value of K we have 
ap = 1.1725(1 ± 5.3%)cH. Analogously to Pioneer 10, this corresponds to H = 76.7 km s~^ 
Mpc~^. Notice that the previous determinations of the initial values for T = 0, corresponds 
with the time when each of Pioneer 10 and 11 begun to follow the hyperbolic orbits. Indeed, 
after Pioneer 10 passed Jupiter and Pioneer 11 passed Saturn, the two spacecrafts followed 
escape hyperbolic orbits (see (Turyshev et al. 2004) and Fig. 2 there). 

Since ^ ~ in case of elliptic orbits, the Pioneer anomaly almost disappears in this 
case. Moreover, ap changes its sign for almost circular orbits. For hyperbolic trajectories ^ 
is much larger. This gives the answer to the question 1 from section 3. 

At about 9.5 AU heliocentric distance ap changes its sign (for K=0) and increases very 
quickly for Pioneer 11, and at about 5.5 AU heliocentric distance (for K=0) ap changes its 
sign and increases very quickly for Pioneer 10 (but this early value of ap is not explored). 
This is the answer to the question 2 from section 3. 

The answer to the question 3 follows also just from the figures 1 and 2. 

Finally, the answer of the question 4 is the following. The formula (3.5) is deduced only 
for the frequencies obtained from the re-transmitted radio signal to the antennae, but not 
for an anomalous acceleration. We saw that the anomalous acceleration studied in section 2 
yields to the following result about the change the distances to the planets. The distances 
to the planets (measured via a laser signal like the distance to the Moon) increase for the 
coefficient 1 + ^tH. This change is not sufficient to be detected for the planets. For the 
distance Earth-Moon it is 1.83 cm per year, while the measured increment is 3.8 cm per 
year. The rest part of 1.97 cm per year is a result of the tidal dissipation. About the 
other measurements in the Earth-Moon system and the explanation according to the time 
dependent gravitational potential see (Trencevski 2006). 
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Fig. 1. — Represents the function Qp/cH for Pioneer 10 as a function of the distance R 
measured in AU. K=15 corresponds that measuring the time started 235.5 days before the 
spacecraft was at the perihehon. 
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Fig. 2. — Represents the function ap/cH for Pioneer 11 as a function of the distance R 
measured in AU. K=30, 45, 60 corresponds that measuring the time started 323, 485, 647 
days before the spacecraft was at the perihehon. 



